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ABSTRACT 

The rotational transitions of carbon monoxide (CO) are the primary means of investigating the density 
and velocity structure of the molecular interstellar medium. Here we study the lowest four rotational 
transitions of CO towards high-latitude translucent molecular clouds (HLCs). We report new obser- 
vations of the J = (4-3), (2-1), and (1-0) transitions of CO towards eight high-latitude clouds. The 
new observations are combined with data from the literature to show that the emission from all ob- 
served CO transitions is linearly correlated. This implies that the excitation conditions which lead to 
emission in these transitions are uniform throughout the clouds. Observed 13 C0/ 12 C0 (1-0) integrated 
intensity ratios are generally much greater than the expected abundance ratio of the two species, indi- 
cating that the regions which emit 12 CO (1-0) radiation are optically thick. We develop a statistical 
method to compare the observed line ratios with models of CO excitation and radiative transfer. This 
enables us to determine the most likely portion of the physical parameter space which is compatible 
with the observations. The model enables us to rule out CO gas temperatures greater than ~ 30 K, 
since the most likely high-temperature configurations are 1 pc-sized structures aligned along the line of 
sight. The most probable solution is a high density and low temperature (HDLT) solution, with volume 
density, n = 10 4 ' 5±0 ' 5 cm -3 , kinetic temperature, Tk ~ 8K, and CO column density per velocity interval 
N co /AV = 10 16 - 6±0 - 3 cur 2 /( kmr 1 ). The CO cell size is L ~ 0.01 pc (~ 2000AU). These cells are 
thus tiny fragments within the ~100 times larger CO-emitting extent of a typical high-latitude cloud. 
We discuss the physical implications of HDLT cells, and we suggest ways to test for their existence. 

Subject headings: ISM: clouds — ISM: molecules — submillimeter 



1. INTRODUCTION 

The rotational transitions of carbon monoxide (CO) are 
the most important probes used to trace the density and 
velocity structure of molecular clouds. Due to the low 
electric dipole moment of the CO molecule, the J = (1-0) 
and (2-1) transitions are easily excited into emission in 
molecular clouds with modest extinctions (A v J> 1 mag- 
nitude). Because of this and the fact that the frequen- 
cies of the (1-0) and (2-1) transitions lie in two relatively 
transparent atmospheric windows, they are the most com- 
monly observed molecular emission lines in the interstellar 
medium (ISM). Unfortunately the information given by 
observing these transitions is limited, since they are not 
sensitive to volume densities above the critical density for 
collisional excitation, n crit ~ 10 3 cm~ 3 . To probe gas at 
higher density demands that transitions with higher ener- 
gies be observed. Star-formation regions are bright in the 
high-J CO transitions, but modelling these transitions is 
difficult because (1) the large CO column density in such 
regions renders their interiors completely opaque, and (2) 
randomly distributed energy sources (stars) increase dra- 
matically the complexity of the problem. Modelling suc- 



cess is more likely in regions which are relatively optically 
thin and far from star formation. 

The simplest class of molecular cloud observable in CO 
is the translucent cloud (van Dishoeck & Black 1988). 
Translucent clouds have modest visual extinctions (A v w 
1 — 5 mag), and therefore are intermediate between low- 
extinction diffuse clouds, which contain few molecules ob- 
servable in emission, and optically thick dark clouds. The 
high Galactic latitude (b ;> 15°) molecular clouds (HLCs; 
see review by Magnani, Hartmann, & Speck 1996, and ref- 
erences therein) are a set of nearby ((d) ~ 105 pc; Magnani 
et al. 1996) translucent clouds associated with IRAS cir- 
rus (Low et al. 1984). The HLCs are generally far from 
sources of far-ultraviolet (FUV) (6eV < hv < 13.6 eV) 
radiation and most are probably exposed to the average 
interstellar FUV radiation field (see van Dishoeck & Black 
1988; Draine 1978; Habing 1968). High-latitude clouds are 
photochemically and energetically the simplest molecular 
objects in the ISM. They have equal amounts of carbon 
gas in atomic form, C I, as in CO, with an average C/CO 
column density ratio of ~ 1 (Stark & van Dishoeck 1994; 
Ingalls, Bania, & Jackson 1994; Ingalls, et al. 1997). They 
are essentially isolated photodissociation regions (Tielens 
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& Hollcnbach 1985) whose structure is dominated by the 
transition from atomic to molecular gas which occurs on 
the surface of all molecular clouds. Most HLCs are uncom- 
plicated by internal energy sources. Some of the HLCs 
contain T Tauri stars (Magnani, Caillault, Buchalter, & 
Bcichman 1995; Pound 1996) and others have marginally 
bound dense cores (Reach et al. 1995), but on the whole 
HLCs are minimally affected by stars. Thus they are the 
most basic clouds observable in CO emission and are ideal 
candidates for multi-transition CO studies. 

The most comprehensive study of CO in translucent and 
high-latitude clouds was done by van Dishoeck, Black, 
Phillips, & Gredel (1991). They observed 4 transitions of 
12 CO and 13 CO in a sample of <~ 25 clouds. They derived 
relatively low gas volume densities (n — 200 — 5000 cm~ 3 ) 
and high kinetic temperatures (Tk > 20 K), but they 
warned that their data are also compatible with higher 
densities and lower temperatures. An ongoing investiga- 
tion of the physics and chemistry of translucent clouds 
(Turner 1994; see also Turner, Terzieva, & Herbst 1999 
and references therein) supports the low density, high tem- 
perature conclusion. On the other hand, a study of CO 
emission towards the translucent edge of a cloud in the 
Perseus-Auriga complex (Falgarone & Phillips 1996, here- 
after FP) claimed that the CO emission is produced by 
beam-diluted clumps with high densities (n > 10 4 cm~ 3 ) 
and low temperatures (T k <; 15 K) (see also Falgarone, 
Phillips, & Walker 1991). Interestingly, FP admitted that 
the observations could also be produced by a more homo- 
geneous CO medium with lower densities (n ~ 200 cm~ 3 ) 
and higher temperatures (Tk ~ 30 — 50 K) . 

In general the solution of the statistical equilibrium 
equations rests on an assumption about either the kinetic 
temperature or the cloud geometry (van Dishoeck et al. 
1991). Clearly such solutions arc not unique, since colli- 
sional excitation rates are affected by both temperature 
and density. Observations of the rotational transitions of 
high dipolc-momcnt molecules in translucent and high 
latitude clouds (Reach, Pound, Wilner & Lee 1995; Hei- 
thausen, Corneilussen, & Grossman 1998) support the as- 
sertion made by FP that significant quantities of dense gas 
exist, but they do not prove that most of the CO emission 
originates in such gas. Even the high-density conclusions 
of Reach et al. (1995), and Heithausen et al. (1998) rely 
on assumed low Tk values (5-20 K). 

The CO emission towards translucent clouds shows some 
surprising regularities which should be exploited in a mul- 
titransition study. For example, the 12 CO (2-1) intensity 
is linearly proportional to the (1-0) intensity in clouds 
without star formation, for all positions and velocities 
observed (FP; Falgarone et al. 1998). This important 
fact implies that the excitation conditions which lead to 
emission in these transitions are uniform throughout the 
clouds. FP argued that the simplest situation which would 
produce uniform excitation conditions was the high den- 
sity, low temperature (HDLT) one (n > 10 4 cm~ 3 ;Tk ~ 
10 K), since the CO level populations would be thermal 
under these conditions. They eliminated the low density, 
high temperature (LDHT) case (n < 10 4 cm~ 3 ; T k £ 20 K) 
because the implied line-of-sight structure sizes (column 
density divided by volume density) are greater than the 
~ 0.06 pc size of the smallest structure observed in pro- 



jection on the plane of the sky. Implicit in their argument 
is the assumption that the size of a structure along the 
line-of-sight equals (on average) the size of the smallest 
observed structure in the plane of the sky. 

This may not be strictly true. The best direct evidence 
for the existence of tiny CO clumps in translucent clouds 
is the population of unresolved <~ 2000 AU (0.01 pc )-sizcd 
structures seen in velocity channel maps of the first IRAM 
key-project (Falgarone et al. 1998). These dense unre- 
solved cores are not isolated spheres, however, but usually 
have the shape of elongated filaments which are unresolved 
in one direction only. Thus there is little reason to sup- 
pose that the average line of sight size of CO emission cells 
equals the smallest observed plane-of-sky structure size. 

Given this geometrical uncertainty, data for many tran- 
sitions are needed to decide if the primary source of CO 
emission is gas with LDHT or HDLT conditions. In partic- 
ular, higher energy transitions of CO should be observed 
towards translucent clouds. Although they observed the 
(4-3) transition of 12 CO towards some positions in the 
cloud they studied, FP were not able to determine if ex- 
citation of this transition was uniform, as it was for the 
(1-0) and (2-1) transitions. If the CO level populations 
are thermal, then the discovery of uniform J = 4 excita- 
tion might support their claim of the ubiquity of gas with 
n > 10 4 cm~ 3 , since the critical density for thermal exci- 
tation of the (4-3) transition is <~ 10 4 ~ 5 cm~ 3 . 

In this paper we report observations of the J = (4-3) ro- 
tational transition of 12 CO towards a sample of 8 translu- 
cent high-latitude molecular clouds at southern declina- 
tions (§2 and §3). Our new data are used in conjunc- 
tion with the 12 CO and 13 CO observations made by van 
Dishoeck et al. (1991) to constrain the physical conditions 
in the HLC regions emitting CO radiation (§4). We show 
that the emission from all CO transitions up to (4-3) is 
linearly correlated, i.e., the excitation is uniform. We use a 
statistical argument to demonstrate that inclusion of our 
new (4-3) data narrows significantly the allowed region 
of parameter space, making the HDLT solution the most 
likely one. We discuss the physical implications of gas with 
HDLT properties, particularly thermal pressure imbalance 
with the ambient interstellar medium (§5). Finally we sug- 
gest a method to test for the existence of HDLT gas. 

2. OBSERVATIONS 
2.1. The source sample 

Eight HLCs at southern declinations were observed (see 
Table 1 for Galactic identifications). Seven of the clouds 
were listed in the catalog of Keto & Myers (1986). The ob- 
served source positions are listed in Ingalls et al. (1997). In 
addition, a previously unknown translucent high latitude 
cloud, G259.4-13.3 [a(1950.0) = 7 /l 29 m 42 s ; (5(1950.0) = 
— 47°00'41"], was observed. This cloud was discovered 
during a search made towards IRAS cirrus clouds for new 
12 CO (2-1) sources (Ingalls 1999). 

The observations described here were part of a study of 
the overall gas-phase carbon content in translucent clouds 
(Ingalls, Bania, & Jackson 1994; Ingalls et al. 1997; Ingalls 
1999). The seven clouds from the Keto & Myers (1986) 
catalog were each recently detected in neutral atomic car- 
bon [C i] emission by Ingalls et al. (1997). Prior to this, 
most of the clouds had only been observed in CO at 8. '7 
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angular resolution (Keto & Myers 1986). Cloud G225.3- 
66.3, has also been observed at higher resolution in CO 
emission and in the photographic B band by Stark (1995). 

Cloud G316. 5+21.0 has a low mass star embedded in it 
(Ingalls et al. 1997). Since this cloud has its own inter- 



nal energy source it may be hotter and more dense than 
the typical HLC. A comparison with the more quiescent 
clouds in the sample may provide insight into the effects 
of low mass star formation on translucent gas. 



TABLE 1 

CO (4-3) and (2-1) Properties of Southern High-Latitude Clouds 

12 CO (J = 4-3) a 12 CO (J = 2-l) b 



Cloud V L 

SR C T MB jT MB dv d 

(kms -1 ) (K) (Kkms- 1 ) (kms -1 ) (K) (Kkms -1 ) 



G225.3-66.3 


0.4 


0.54±0.15 


0.43 ±0.37 


0.6 


1.4±0.2 


2.8 ±0.4 


G295.3-36.2 


2.1 


0.64±0.13 


0.62±0.12 


1.8 


3.0 ±0.5 


6.6±1.0 


G259.5-16.5 


3.6 


1.41 ±0.22 


1.3 ±0.3 


3.5 


5.6 ±0.8 


5.7±0.9 


G259.4-13.3 


5.0 


2.2 ±0.3 


2.6 ±0.4 


4.9 


8.1 ±1.2 


22.8 ±3.4 


G292.7-19.8 


4.6 


0.34 ±0.06 


0.2 ±0.2 


4.4 


3.7 ±0.6 


7.9 ±1.2 


G294.4-14.3 


4.5 


0.31 ±0.07 


0.64 ±0.12 


4.3 


3.8 ±0.6 


8.0±1.2 


G272.9+29.3 




< 0.2 e 


< 0.24 e 


-3.2 


2.1±0.3 


2.4 ±0.4 


G316.5+21.0 


-6.2 


3.2 ±0.5 


6.4 ±1.0 


-6.5 


13.6±2.0 


33. ±5. 



a Based on 12 CO (4-3) map data averaged over a 5' diameter circle. 
b Single pointing at the nominal cloud position with the AST/RO 5' beam. 

c Determined from a Gaussian fit. A 15% calibration uncertainty is assumed for T MB ; this is added in quadrature 
with AT RMS to give the quoted uncertainty. 

integrated intensity. Errors are determined by adding in quadrature the measurement uncertainty (based on 
RMS) and 15% calibration uncertainty. 

e Upper limits are 3xAT RMS for T MB and 3xAT RMS x 1.8 kms -1 for jT ME dV. 



G225.3-66.3 




Fig. 1. — Observations of high-latitude clouds. The greyscale maps show HIRES-processed IRAS lOO^tm emission. Black contours represent 
AST/RO 12 CO (J = 2-1) integrated intensity. A 5' circle is drawn on each map in white, at the nominal source position, defined by Keto & 
Myers (1986). The average 12 CO (2-1), 12 CO (4-3), and [C I] ( 3 Pi -> 3 P ) spectra measured in this circle with AST/RO are shown on the 
right. Gaussian fits are superimposed on the spectra. The C I data were originally published by Ingalls et al. (1997). The 12 CO (4-3) and C I 
intensities have been multiplied by 3. (a) Cloud G225.3-66.3. This source was not mapped in 12 CO (2-1) emission. (6) Cloud G295.3-36.2. 
Contours of 12 CO (2-1) emission range from 2.2 to 2.8 K kms -1 , in 2.2 K kms -1 intervals. 
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FlG. 2. — Observations of high-latitude clouds. See caption to Figure Mfor description, (a) Cloud G259.5-16.5. Contours of 12 CO (J = 
2—1) emission range from 1.9 to 7.6 K kms -1 , in 1.9 K kms -1 intervals, (b) Cloud G259.4-13.3, discovered by Ingalls (1999), with nominal 
source position taken from that work. Contours of 12 CO (2—1) emission range from 7.7 to 30.8 K kms -1 , in 7.7 K kms -1 intervals. 



2.2. AST/RO Observations of CO (4~3) and (2-1) 
emission 

Observations of the 461.041 GHz (J = 4-3) transition of 
12 CO were made during the period 1997 June 29 through 
July 13, using the Antarctic Submillimctcr Telescope and 
Remote Observatory (AST/RO). AST/RO is a 1.7m diam- 
eter telescope located at the National Science Foundation 
Amundsen-Scott South Pole Station (Stark et al. 1997). 
The telescope efficiency at 461 GHz was estimated from 



skydip measurements (Chamberlin, Lane, & Stark 1997) 
to be r\i — 0.77. As mentioned by Ingalls et al. (1997), 
r\i is close to the main-beam efficiency, ij mb . The 461 
GHz transition was measured using the AST/RO 492/810 
SIS waveguide receiver which had a ~ 130 K double side- 
band receiver noise temperature. The intermediate fre- 
quency output was sampled by a 2048 channel, 1.4 GHz 
bandwidth acousto-optical spectrometer (AOS) (Schieder, 
Tolls, & Winnewisser 1989). The velocity resolution was 
^V 46 i = 0.62 kms" 1 channel -1 . System noise tempera- 




FlG. 3. — Observations of high— latitude clouds. See caption to Figure |l| for description, (a) Cloud G292.7— 19.8. Contours of 12 CO (J = 
2-1) emission range from 2.8 to 9.9 K kms -1 , in 2.4 K kms -1 intervals, (b) Cloud G294.4-14.3. Contours of 12 CO (2—1) emission range 
from 2.4 to 9.9 K kms -1 , in 2.5 K kms -1 intervals. 
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Fig. 4. — Observations of high-latitude clouds. The data shown here correspond to those described in the caption to Figure hj, with the 
addition of FCRAO 12 CO (1-0) data in part a. (a) Cloud G272.9+29.3. White contours of 12 CO (J = 1-0) emission range from 1.9 to 
3.9 K kms -1 in 1.9 K kms -1 intervals. The average (1-0) spectrum is superimposed on the plot at right. Black contours of 12 CO (J = 
2-1) emission range from 0.6 to 3.0 K kms _1 in 0.6 K kms" 1 intervals. The 12 CO (4—3) line was not detected from this source and the 
nondetection spectrum is shown. (6) Cloud G316. 5+21.0. Contours of 12 CO (2-1) emission range from 13.8 to 34.5 K kms -1 , in 6.9 K kms -1 
intervals. 



tures were of order 1800-3000 K . The half power beam 
width of the AST/RO telescope at 461 GHz was approxi- 
mately 180". 

Observations of the 230.538 GHz (J = 2-1) transition of 
12 CO were made in 1996 November using the AST/RO 230 
GHz SIS receiver which had a double sideband noise tem- 
perature of ~ 500 K. Based on measurements of r\i from 
skydips (see above), the value r] MB = 0.54 was adopted for 
230 GHz. The same AOS used for the (4-3) measurements 
was used for the (2-1) data, yielding a velocity resolution 
of <5V 230 = 1.34 kms -1 channel -1 . The system noise tem- 
perature at 230 GHz was ~ 1800 K. The beamsize of the 
AST/RO telescope at 230 GHz is - 5'. 

The 12 CO (4-3) and (2-1) observations were obtained 
by position-switching 1° in azimuth. Because of the lo- 
cation of AST/RO, position-switched offsets were both at 
constant elevation and at fixed positions on the sky. Inten- 
sity calibration was accomplished by measuring two black- 
body loads at known temperatures (Stark et al. 1997). 
In order to correct for atmospheric attenuation the sky 
brightness temperature was measured at the position of a 
source every ~ 15 minutes. The overall intensity calibra- 
tion was checked by observing periodically the compact 
H II region G291.2-0.8. This procedure yielded intensities 
repeatable to within ±15% at both 230 GHz and 461 GHz. 

2.3. FCRAO observations of CO (1-0) towards cloud 
G272.9+29.3 

The HLC G272.9+29.3 was observed with the Five Col- 
lege Radio Astronomy Observatory (FCRAO), in Amherst 
MA, on 1996 February 17 through 18. The FCRAO 14m 
telescope was used to map the 12 CO (1-0) transition at 
115.2712 GHz. The telescope has a beamsize of 45" and 
a 115 GHz main beam efficiency of ?7mb 

« 0.36 at 20° 



elevation (Ladd & Heyer 1996). The observations were 
carried out in frequency-switched mode, with the signal 
and reference spectra spaced 4.5 MHz apart. The spectra 
were measured with the QUARRY focal plane array re- 
ceiver (Erickson et al. 1992), and the data were sampled 
by an array of 15 autocorrelation spectrometers with 1024 
channels, separated by 77 kHz per channel, giving a veloc- 
ity resolution of <5V 115 = 0.2 kms -1 channel -1 . Since the 
source is a Southern hemisphere object with declination 
5 = —26°, its elevation was rather low, ranging from 16° 
to 19°. The single-sideband system noise temperature was 
approximately 2000 K for these observations. 

3. OBSERVATIONAL RESULTS 
3.1. CO (4-3) and (2-1) Properties of HLCs 

Observed 12 CO (4-3) and (2-1) spectra are plotted 
alongside maps of each source in Figures [j] to [|. The (4- 
3) transition was detected towards 7 of 8 clouds observed. 
The 12 CO (4-3) HLC data consist of half-beam sampled 
measurements made ona~ 5' x 5' grid, in a manner similar 
to the C I "superbeam" observations described by Ingalls 
et al. (1997). These maps were averaged over the 5'-square 
grids to produce the spectra shown in Figures @to|. The 
C I "superbeam" average spectra taken towards the same 
positions (Ingalls et al. 1997) are also displayed in the fig- 
ures. Intensities are main-beam brightness temperatures, 
T MB , i.e., Rayleigh- Jeans antenna temperature, Ta, cor- 
rected for both atmospheric attenuation and main-beam 
efficiency (T MB = TX/i] MB ). 

Each source was observed in 12 CO (2-1) emission at its 
nominal center with the ~ 5' AST/RO beam. For seven of 
the sources, these spectra are part of ~ 30' x 30' half-beam 
sampled maps, shown as black contours overlaid on IRAS 
100 /im greyscale images in Figures to |. The IRAS im- 
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ages were all processed using the HIRES technique (Levine 
et al. 1993), and have angular resolution ~ 100". The cali- 
bration uncertainty for the HIRES images is of order 20%. 

Gaussian fits were made to the CO (4-3) and (2-1) spec- 
tra and are shown superimposed on these spectra in Fig- 
ures [| to ||. Table 1 lists the line center velocity, V LSR , and 
the peak intensity, T MB , from the Gaussian fits. Also listed 
is the integrated line intensity, J T MB dV . Error bars for 
both transitions were determined by adding in quadrature 
the instrumental noise (AT RMS ) and a calibration uncer- 
tainty of 15%: 



ct(T mb ) = y (AT RMS ) + (0.15T MB )' 



(1) 



Gaussian fit linewidths, AT^, range from 0.9 to 3 kins" 1 
with an average of (1.7 ±0.5) kms -1 , slightly greater than 
the velocity resolution of the (2-1) observations, so the 
lines are only marginally resolved. For this reason both 
T MB and AV are poorly-determined. We will therefore fo- 
cus our attention on the integrated intensity, J T MB dV, 
which is a better defined quantity. Since (4-3) radia- 
tion was not detected towards G272.9+29.3, we estimate 
a 3cr upper limit to the (4-3) integrated intensity for this 
source using the average linewidth: J* T MB (4 - 3) dV < 
3 x AT R ms x 1.7 kms -1 . 

3.2. The CO (1-0) Map of G272.9+29.3 

The (1-0) transition of 12 CO was mapped in source 
G272.9+29.3 at 45" angular resolution. The map con- 
sisted of a 29 x 11 grid of pixels, evenly spaced at ~ 25", 
one-half the FCRAO beamwidth. The resulting map of 
12 CO (1-0) integrated intensity is superimposed as white 
contours on Figure |^ (a). The average spectrum in a 5' 
"superbeam" at the nominal source center is plotted on 
the right side of this figure. This spectrum was resolved in 
velocity by the FCRAO spectrometer: a Gaussian fit gives 
the line center velocity V LSIi = —3.1, linewidth AV = 1.0 
and peak intensity T MB = 5.6 ± 1.8. 

3.3. The HLCs as Representative CO Clouds 

The HLCs are the average Galactic high-latitude 
sources of CO emission. This can be shown in two ways. 
First, we can use the all-sky CO survey done with the Far 
Infrared Absolute Spectrometer (FIRAS) instrument on 
board the Cosmic Background Explorer (COBE) satellite. 
This instrument was used to observe nearly the entire sky 
at 7° resolution in the first six transitions of 12 CO (Wright 
et al. 1991; Bennett et al. 1994). The data are available 
on the World Wide Web. 3 The average value of the COBE 
FIRAS high-latitude (|6| > 15°) (4-3)/(2-l) integrated in- 
tensity ratio is ((4 - 3)/(2 - 1)) CO be = 0.84 ± 0.06. The 
intensities were reported in units of erg cm -2 s - 1 sr~ 1 . The 
FIRAS integrated intensity ratio may be compared with 
the ratio of integrated Rayleigh- Jeans brightness temper- 
atures measured here for HLCs if one uses the conversion 
formula: 



The function Q(u) = 2 x I0 5 kis 3 /c 3 : Q(230) = 1.25 x 10" 8 
and Q(461) = 10~ 7 ergcm- 2 s" 1 sr -1 /(K kms" 1 ]. The 
average HLC value measured here is [Equation p[ mul- 
tiplicd by Q(461)/Q(230) = 8] ((4 - 3)/(2 - 1)J HLC = 
0.88±0.16, which is nearly identical to the COBE average. 

The infrared properties of the HLCs studied here can 
also be compared with the average properties of interstel- 
lar cirrus clouds. HIRES-processed images of both 60 and 
100/im emission observed with IRAS were obtained for 
each cloud in the sample. The 60/Ltm surface brightness, 
I 60 , is plotted as a function of the 100/zm surface bright- 
ness, I 100 , in Figure |. Both the 60 and lOO^rn measure- 
ments are averages over the 5' "superbeam" centered on 
each source, described in §3.1. The Diffuse Infrared Back- 
ground Experiment (DIRBE) on board the COBE satel- 
lite was used to observe the infrared spectrum of the entire 
Galaxy. The average result for all Milky Way gas corre- 
lated with H I emission, (Lo/Loo)cobe — 0.160, was derived 
by Dwek et al. (1997). This is the same as the mean value 
for the HLCs studied here, (WLoo)hlc = 0.156 ± 0.016 
(the authors did not quote an error for this ratio). A 
straight line indicating I 60 = 0.16I loo is superimposed on 
Figure S. Notice that most of the clouds fall on this 
line, with the exception of source G316. 5+21.0 which has 
Ieo = 77+15 MJysr" 1 . For this cloud I 60 /I 100 = 0.44±0.12, 
which is much larger than that of the other clouds in the 
sample. This implies that G316. 5+21.0 has a greater rel- 
ative abundance of the "very small grains" which produce 
most of the 60/im emission (e.g., Desert, Boulanger, & 
Puget 1990). The low mass star which is embedded in 
the source must be responsible for this, presumably by 
heating the dust in excess of the interstellar radiation field 
(ISRF) and evaporating the mantles of larger grains which 
produce 100 /xm emission. 

3.4. Correlations Between Integrated Intensities 

Here we give evidence for uniform CO excitation in 
translucent clouds. Figure ^ shows the observed (4-3) and 
(2-1) integrated intensities for the HLC sample, plotted as 
a function of I 100 . Both graphs in Figure ^ show the same 
behavior, although the trend is not necessarily monotonic. 
Figure is a plot of /(4-3) = / T MB (4 — 3) dV as a func- 
tion of J (2 _ 1) = /T MB (2 - l)dV. The (4-3) integrated 
intensity is linearly correlated with the (2-1) integrated 
intensity. The correspondence is significant, with a cor- 
relation coefficient of r = 0.95. Linear regression to the 
data, weighted by the errors, gives: 



7 (4-3) 



(0.11 ± 0.02) J ( 2_i 



(3) 



A line representing Equation || is superimposed on the 
data in Figure fj]. The 1(4,-3) -intercept was held fixed at 
zero, i.e., we have assumed that there is no (4-3) emission 
when there is no (2-1) emission, and vice-versa. Source 
G316. 5+21.0 has an anomalously high (4-3)/(2-l) inte- 
grated intensity ratio. As explained in §3.3, the star em- 
bedded in this cloud is likely to make the gas hotter than in 
the average translucent cloud, whence high level CO tran- 
sitions will be excited more readily into emission. Since 
the linear fit to the (4-3) and (2-1) data was weighted by 

3 The Web site is http://www.gsfc. nasa.gov/astro/cobe/ . The COBE datasets were developed by NASA's Goddard Space Flight Center 
under the guidance of the COBE Science Working Group and were provided by the National Space Science Data Center. 
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Fig. 6.— Observed 12 CO (J = 4-3) {Top) and (J = 2-1) (Bottom) properties of HLCs. The integrated main beam brightness temperature, 
J TmboIV, is plotted here for each transition as a function of IRAS 100/im surface brightness. 



the errors, this source (which had rather large error bars) 
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did not contribute significantly to the fit. 

The linear relationship between 1(4-3) and ^(2-1) (Fig- 
ure pi) is reminiscent of the correlation between the 12 CO 
(2-1) and (1-0) brightness temperatures seen towards 
"pre-star-forming regions" by Falgarone et al. (1998). 
There is nothing special about the HLC positions we 
observed; they are the nominal "cloud center" positions 
(Kcto & Myers 1986; Ingalls 1999), and do not necessar- 
ily coincide with the peak of CO emission. Recall that 
all quantities quoted in Table 1 refer to measurements in 
a 5' circle (see Figures |l| to ^) made towards these cen- 
ters. These nominal positions have 100/im surface bright- 
nesses ranging over 2 orders of magnitude, so the sample 
is obviously not homogeneous. Nevertheless, the ratio be- 
tween the (4-3) and (2-1) integrated intensities, as well 
as between the 60 and 100/im intensities, is remarkably 
uniform. The simplest interpretation is that the CO emis- 
sion from all translucent clouds originates from gas with 
the same physical conditions, even though different lines 
of sight contain different amounts of such gas. This was 
the conclusion of Falgarone et al. (1991) and FP, based on 
their study of the edge of a cloud in the Perseus complex. 

A linear correlation seems to exist between all pairs of 
CO transitions observed towards HLCs. The most exten- 
sive database of CO observations of HLCs is that of van 
Dishoeck et al. (1991). Figure || shows -f(i-o) (Top) and 
1(3-2) (Bottom) plotted as a function of f° r the ~ 25 

translucent and/or high-latitude positions detected by van 
Dishoeck et al. (1991; their Table 7). Both the I(i-o) vs - 
/(2-1) an d the 1(3-2) vs - ^(2-1) datasets have a highly sig- 
nificant correlation (r = 0.94 for each dataset), of which 



van Dishoeck et al. (1991) were apparently unaware. Lin- 
ear regression to the data gives the following fits: 

7 (1 _ 0) = (1.3±0.4)7 (2 _ 1); (4) 

7(3-2) = (0.55 ± 0.16) 7(2.!). (5) 

As for Equation ||, the fits were constrained to have an 
intercept value of zero. We weighted the data using 30% 
error bars, corresponding to the overall calibration uncer- 
tainty quoted by van Dishoeck et al. (1991). This cor- 
relation also applies to the cloud sample used by Turner 
(1993). For Turner's translucent cloud sample, the aver- 
age 12 CO (l-0)/(2-l) integrated intensity ratio for nine 
sources is 1.6 ± 0.3, which is indistinguishable (within the 
errors) from the result given above. 

The slope in Equation |] is equivalent to the constant 
(2-l)/(l-0) brightness temperature ratio observed by Fal- 
garone et al. (1998). Considering each spectral channel 
individually, they measured T( 2 -i)/T(i-o) = 0.65 ± 0.15 
for 80% of the data points in their maps made towards 
three translucent clouds. In comparison, the reciprocal of 
the Equation | slope is 7 (2 _ 1) /7( 1 _ ) = 0.77 ± 0.24. Note 
that Falgarone et al. (1998) and FP did not study iso- 
lated translucent objects, but since they obtained similar 
values for I(2-i)/ I(i-o) it is possible that the linear rela- 
tionship is also a characteristic of the translucent edges of 
non-starforming dark clouds. 

This result is not limited to translucent regions. A large 
scale survey of the Milky Way molecular ring has obtained 
an average value of the (2-l)/(l-0) ratio of 0.85 ± 0.63 
with a median of 0.69 (Chiar, Kutner, Verter, & Leous 




Fig. 7. — Integrated main beam brightness temperature Ico (4-3) versus Ico (2-1)- The straight line fit, ^(4-3) = 0.11/(2-1) 1 i s 
superimposed. 
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Fig. 8. — Integrated intensity of the 12 CO (1-0) (Top) and (3-2) (Bottom) transitions, plotted as a function of the (2-1) integrated intensity 
for translucent and high-latitude clouds (from Table 7 of van Dishoeck et al. 1991). The linear fits to these data, J(i_o) = 1-3/(2-1) an d 
^(3-2) = 0.55/(2_i), are superimposed. 



1994). Similarly, the mean value for the entire Galac- 
tic plane, (T( 2 _i)/T( 1 _ )) = 0.77, was derived by Sanders 
et al. (1993; also see Sanders et al. 1986), and Handa et al. 
(1993) measured a mean ratio of 0.7 for the Galactic molec- 



ular ring (also see Dame et al. 1986). Thus essentially all 
surveys of CO consistently show a linear correlation be- 
tween (2-1) and (1-0) integrated emission. 

This correlation also exists for the CO peak brightness 



6 F 




Fig. 9. — Integrated 13 CO (1-0) intensity versus integrated 12 CO (1-0) intensity for translucent and high-latitude clouds (Table 7 of van 
Dishoeck et al. 1991). The linear fit 7( 13 CO) = 0.13/( 12 CO) is superimposed as a solid line. The upper limit behavior expected for optically 
thin emitting regions, J( 13 CO) = 0.028/( 12 CO) (see text), is superimposed as a dashed line. In general the 12 CO (1-0) emission is optically 
thick. 
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temperature. Peak brightness temperatures were listed by 
van Dishoeck et al. (1991) for the translucent sources with 
spectrally resolved (1-0), (2-1), and (3-2) line profiles. We 
find that these are also correlated, and linear fits give the 
same slopes as the integrated intensity fits. In general, for 
any two of these transitions, "1" and "2" , 



s(l)dV 
3 {2)dV 



Tmb(1) 

Tmb (2) 



(6) 



This is a trivial statement for Gaussian line profiles but 
according to van Dishoeck et al. (1991) not all of their 
spectra were well-described by a Gaussian shape. There- 
fore, linear correlation of integrated intensities is not solely 
due to a correlation in linewidths which might occur in op- 
tically thick gravitationally bound systems with different 
masses. 

We make one last comparison, that between the inte- 
grated intensities of the (1-0) transition of the 13 CO and 
12 CO isotopomers. The data are listed in Table 7 of van 
Dishoeck et al. (1991) and plotted here in Figure |j| For 
these data the correlation coefficient, r — 0.64, so there is 
only an 80% chance that the data are correlated (Taylor 
1997). Linear regression (with zero intercept) yields: 



J( 13 CO) 



(0.13±0.04)/( 12 CO). 



(7) 



The fit is not very robust, but the slope 0.13±0.03 gives 
a good estimate of the average 13 CO/ 12 CO intensity ra- 
tio. (We obtain a comparable value for this ratio using 
the 13 CO and 12 CO HLC data of Keto & Myers 1986). 
This is approximately five times higher than what is ex- 
pected if the CO-emitting regions are optically thin to both 
13 CO and 12 CO radiation, assuming the abundance ratio 
[ 12 C]/[ 13 C] = 60 (Boreiko & Betz 1996; Langer & Penzias 
1993). Taking into account isotope-selective photodissoci- 
ation and chemical fractionation, the slope from optically 
thin gas should range from 7( 12 CO)//( 13 CO) w 0.014 to 
0.028 (van Dishoeck & Black 1988). Nearly all of the data 
points fall above the line 7( 13 CO) = 0.028/( 12 CO) which 
is superimposed on Figure |^, so we conclude that most 
of the CO-emitting regions in HLCs are optically thick to 
12 CO (1-0) radiation. 

We have shown that the CO integrated intensities I(i-o) , 
1(2-1) i I(3-2)i an d -^(4-3) are related linearly for a large 
sample of translucent high-latitude molecular clouds. Un- 
fortunately the cloud sample used by van Dishoeck et al. 
(1991) does not overlap with our own sample, so it is un- 
certain whether all of the clouds in both samples fulfil all 
of the linear relationships. This is probably true, however, 
since both samples probe the same kind of gas. By defini- 
tion, all of the clouds are translucent (i.e., A v = 1 — 5, no 
internal UV sources; see Ingalls et al. 1997). In addition, 
they have similar 12 CO (2-1) and 100/xm emission prop- 
erties. For example, the weighted mean 7(2-i)/Iiao ratio 
for our sample is 0.55±0.06K kms -1 ( MJy sr -1 ) -1 . For 
the cloud towards HD 210121, towards which van Dishoeck 
et al. (1991) observed six positions, the average /(2-i)/Iioo 
ratio is 0.6 ± 0.4K kms -1 ( MJy sr -1 ) -1 . It is therefore 
likely that the clouds in both samples are equivalent. In- 
deed, Falgarone et al. (1998) cite the linear relationship 
between intensities of CO rotational transitions as a gen- 
eral characteristic of "clouds of low average column den- 
sity" , i.e., translucent clouds (see also Clemens & Barvainis 



1988; Falgarone et al. 1991; Falgarone, Puget, & Perault 
1992; FP). 

In what follows we interpret the slopes of the linear fits 
(Equations ||, |], [|, an d 0) to be the true (constant) val- 
ues of the respective integrated intensity ratios in the cells 
which emit CO radiation. We take the CO-emitting cells 
to be nearly identical structures which are separated from 
each other in position and/or radial velocity. In this model 
an increase in the integrated intensity of a CO transition is 
attributed solely to an increase in the number of cells being 
observed, i.e., an increase in 0, the beam filling fraction of 
the CO radiation emitting from a particular direction at 
a specific velocity. 



4. MODELLING THE CO-EMITTING REGIONS: A 
STATISTICAL APPROACH 

Here we introduce a statistical method for modelling the 
physical conditions in the CO-emitting cells of translucent 
molecular clouds. As stated above, integrated intensity 
ratios among the first four transitions of CO are constant. 
Define the following integrated intensity ratios for translu- 
cent clouds: 





-1) 


ki- 


-0) 




-2) 




-1) 


'(4- 


-3) 



1(2-1) 



(.Ri±o-i) = 0.77 ±0.24; (8) 



= (R 2 ±a 2 ) = 0.55 ±0.16; (9) 



7[ 13 CO(1-0)] 
7[ 12 CO(l -0)] 



(i?3±cr 3 ) = 0.11 ±0.02; (10) 



(i? 4 ±cr 4 ) = 0.13 ±0.04. (11) 



These ratios are taken from Equations |3[ ^, H, and [?]. Sup- 
pose that fluctuations in the ith ratio, fj, are distributed 
normally with mean Ri and standard deviation a%. The 
probability of measuring the ith ratio to be between and 
n + dri is then: 



p(ri,Ri,ai)dri 



1 



<Tj V 2tt 



<r i -R i r/(2^) dr ._ (12) 



Since the ratios are constant (excluding fluctuations), they 
are statistically independent of each other. The joint prob- 
ability of observing a given set of ratios is thus given by 
the product: 



P(ri,r 2 ,r 3 ,r4)dridr 2 dr 3 dr4 = P(r)dr = [jpfr,,]?,,^)*, 

i=l 

(13) 

This analysis requires that all the CO emission origi- 
nates in gas with the same physical conditions. We of 
course have argued that this is the case (see §3.4). There- 
fore, given a physical model and an input grid of physical 
parameters, the line intensities, T MB , and ratios among 
them, r, can be computed. Equation [Tj}] gives the prob- 
ability distribution in the vector of ratios, r. This can 
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Fig. 10. — Probability analysis of HLCs in (Nco/AV)-temperature space, under LTE conditions (see Table 2). The dashed curve encloses 
the la confidence region when only HLC data for the (3-2) and lower transitions of 12 CO and 13 CO are taken into account. The solid contours 
enclose the 1,2, and 3<r confidence intervals when the new (4-3)/(2-l) data (this paper) are included. The inner, lcr, contour is shaded. The 
units of Nco/AV are cm -2 ( kms" 1 ) -1 . 



be used to estimate the most likely portion of the input 
parameter space which gives the observed ratios. Below 
we evaluate the high density, low temperature (HDLT) 
and the low density, high temperature (LDHT) translu- 
cent cloud scenarios using two simple cloud models, the 
local thermodynamic equilibrium model and the large ve- 
locity gradient model. 

4.0.1. The LTE model 

In the HDLT scenario advocated by FP the CO emis- 
sion originates in cold, high-density regions where the CO 
level populations are in local thermodynamic equilibrium 
(LTE) . This scenario is worth investigating because of its 
extreme simplicity — in LTE the line ratios depend on only 
two parameters, the excitation temperature, T ex , and the 
CO column density per velocity interval, N co /AV. 

The peak intensity of a CO spectral line, T MB , is given 
in LTE by: 



T, 



.53 K 



1 



e (5.53K)J/T ox _ l e (5.53K)J/T BG _ 1 



(14) 



where T BG = 2.7 K is the cosmic microwave background 
temperature and r is the line center opacity. The tran- 
sition is from rotational level J to J — 1. The function 
<f> contains the details of the telescope beam average; if 
the emission originates from gas with uniform properties, 
then 4> is simply the beam filling fraction of the gas. In 



this case, which we assume to hold for HLCs, <j) cancels 
when line ratios are taken. 

We will use the LTE model to predict T MB ratios, and 
then compare them directly with the observed J T MB dV 
ratios (Equation ||) by assuming Gaussian line profiles. 
The LTE opacity for a transition from level J to level 
J — 1 is given by: 



TjJ-l 



8tt> 



3, ,2 



-(2.765 K)J(J+l)/T c; 



3hZ(T cx ) 

r„(5.53K)J/T el! 



1] [N co /AV] 



1.565 x 1CT 15 [ e ( 5 - 53K ) J / T » - i] 



Z(T CX ) 
Nco/AV 



,(2.765 K).7(.7+1)/T c 



kms 



(15) 



Here /i = 0.112 debye is the electric dipole moment of the 
CO molecule and Z(T CX ) is the CO rotational partition 
function: 



Z(T CX ) = £(2J+1) exp[-(2.765K)J(J + l)/T cx ]. 



(16) 

ox , we have corn- 



Using a grid of values of N co /AV and T, 
puted the vector r of ratios under LTE conditions. The 
series in Equation ^ was truncated to 20 terms, since new 
terms added less than 1% to the value of Z. The 13 CO 
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Fig. 11. — Statistical equilibrium calculations for CO lines in clouds with a large velocity gradient (LVG). Each plot maps out observed 
cloud CO line ratios as a function of column density per unit velocity, Nco/AV, and gas volume density for a different temperature. The 
shaded region indicates the range of values of the 12 CO (4-3)/(2-l) ratio reported here for HLCs. Straight lines indicate the total "size" of 
the emitting regions, given by the ratio of column density to volume density, for CO-emitting cells with velocity width AV = 1 kms -1 . The 
units of Nco/AV are cm~ 2 ( kms -1 ) -1 . Density is in units of cm -3 . 




Iog 10 [N C0 /AV] 

Fig. 12. — The probability density field for HLC CO-emitting cells in the LVG model. Slices in the (Nco/AV')— Density plane were made 
for 10 values of the temperature, and are represented as greyscale images. The dashed curves in each plot represent slices through the 3cr 
confidence surface for the sample region of [T^, Nco/AV, n]-space when only HLC data for the (3—2) and lower transitions of 12 CO and 
13 CO are taken into account. The solid contours represent slices through the 3cr confidence surface when the new (4-3)/(2-l) data (this 
paper) are included. The crosses indicate the point in each slice where the probability density xeaches its maximum value, listed as "Peak" 
on each graph. Straight lines indicate the "size" of the emitting regions, as described in Figure 11. 
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intensities were determined in the same way as the 12 CO 
intensities, except the 13 CO abundance was taken to be 
1/60 that of 12 CO. 

TABLE 2 
Predicted CO Line Ratios in LTE a 





(2-l)/(l-0) 


(3-2)/(2-l) 


(4-3)/(2-l) 


5K 


0.61 


0.36 


0.01 


7.5 K 


0.73 


0.59 


0.10 


10 K 


0.79 


0.70 


0.25 


15K 


0.87 


0.79 


0.45 


20 K 


0.94 


0.83 


0.56 


Observed: 


0.77 ±0.24 


0.55 ±0.16 


0.11 ±0.02 



a Column density is assumed to be 5X10 1 



Predicted 12 CO LTE line ratios are listed in Table 2 as a 
function of excitation temperature for N co = 5x 10 16 cm~ 2 
and AV = 1 kms -1 . The bottom row of the table gives 
the actual values observed towards HLCs. The measured 
(2-l)/(l-0), (3~2)/(2-l), and (4-3)/(2-l) ratios are all 
compatible with low temperatures, T cx ps 5 — 10 K. Us- 
ing the probability distribution function, P(r)dr (Equa- 
tion p^[ ), we can determine quantitatively a probability 
distribution in (N co /AV, T ox )-space. Figure [l(] shows 
this. The dashed contour on the figure represents the la- 
confidence interval (68.3% of the total probability) in the 
(Nco/AV, T cx )-plane when only ratios 1, 2, and 4 are 
considered, i.e., when we neglect the new (4-3) data (this 
paper). The solid shaded contour represents the la con- 
fidence interval when the new data are included, and the 
two surrounding solid contours represent the 2 and 3tr con- 
fidence intervals (95.4% and 99.7%, respectively, of the to- 
tal probability). The new data have decreased the size 
of the lcr contour by a factor of ~ 3. The most prob- 
able set of physical properties is: T ox = (8.5 ± 1.5) K 



and N C o/AV = 10 



16.5±0.2 , 



( kms 



Evidently, 



if CO-emitting regions arc in LTE then the temperature 
is low, and the HDLT regime is the only solution. The 
LTE line center opacity (Equation 15) is r = 9.6, 7.2, 2.1, 
and 0.3, for the 12 CO (1-0), (2-1), (3-2), and (4-3) tran- 
sitions respectively. As we deduced from the observations 
(see Figure ||), the CO-emitting regions are optically thick 
to 12 CO (1-0) radiation. In fact, only the (4-3) and higher 
transitions of 12 CO are optically thin. 

4.0.2. The LVG model 

The LTE analysis is necessary but not sufficient to con- 
firm the HDLT scenario because the LTE model does not 
consider the possibility of low density gas. We have sim- 
ply shown that if the gas is dense enough such that the 
transitions are thermalized, then the temperature must be 
low. There is no reason a priori to believe this is true. If 
the temperature is high enough, intermediate levels such 
as J = 4 may be significantly populated even in low den- 
sity gas. Furthermore, it may be radiative trapping which 
causes the thermalization of level populations, rather than 
high densities. The LTE model cannot distinguish between 



such possibilities. 

In this section we calculate the CO level populations 
in statistical equilibrium using the large velocity gradient 
(LVG) formalism (Sobolev 1960; Castor 1971). This ap- 
proximation states that the mean free path of spectral line 
photons, i mfp , is on average much less than the product 
AV (dV/dR)' 1 , where AV is the total cloud velocity dis- 
persion and dV I dR is the cloud velocity gradient which is 
a systematic function of position. In the LVG approxima- 
tion different regions in a cloud are effectively decoupled 
in velocity space, whence line radiation emitted from one 
region of the cloud does not interact with any other re- 
gion. This simplifies considerably the radiative transfer 
computation. 

Of course, the LVG approximation does not strictly hold 
in HLCs — it is unlikely that the velocity field in HLCs is 
so well-ordered as to produce a systematic velocity gradi- 
ent. In general the clouds are not gravitationally bound 
(MBM; Keto & Myers 1986; Pound & Blitz 1993; Hei- 
thausen 1996), nor do they possess evidence of large-scale 
outflows or rotation (e.g., Pound et al. 1991). Instead, 
the velocity field in these clouds is most likely turbulent 
(Falgarone 1994; Falgarone et al. 1994), hence the veloc- 
ity is a random function of position. Nevertheless, Z mfp is 
probably much less than the velocity correlation length of 
the turbulence, i.e., the spectral lines form under macro- 
turbulent conditions [see for example the simulations of 
Park & Hong (1995); see also the review in Falgarone 
et al. (1998)]. Macro-turbulence prevents distinct emis- 
sion regions of a cloud from "shadowing" each other in 
velocity space, which is statistically equivalent to a large 
systematic velocity gradient in its effect on the radiative 
transfer through the cloud. 

Since the level populations both determine and are de- 
termined by the radiation field, the solution of the equa- 
tion of radiative transfer for CO line radiation with self- 
consistent level populations was performed numerically. 
The level populations in a statistical steady state, and the 
corresponding radiation field were computed for a grid of 
clouds, each with constant volume density, n, kinetic tem- 
perature, Tk, and CO abundance, X co - At the surface 
the incident radiation field was taken to be that due to 
the 2.7 K cosmic microwave background. It was further- 
more assumed that there is no internal radiation source 
in the clouds. Collisional rate coefficients for CO exci- 
tation by H2 were those computed by Flower & Lau- 
nay (1985). A constant abundance X co = 8 x 10 -5 
was assumed, i.e., in the CO-emitting regions approxi- 
mately 80% of the carbon gas is in CO. This is based 
on the local ISM gas-phase carbon abundance which is 
observed to be [C]= (1.4 ± 0.2) x 10~ 4 (Sofia, Cardelli, 
Gucrin, & Meyer 1997). The results do not change dras- 
tically if X co varies between 5 and 10 xl0~ 5 . The abun- 
dance of 13 CO was assumed to be 1/60 that of 12 CO, i.e., 
[ 13 CO]/[ 12 CO] = [ 13 C]/[ 12 C]. We have ignored the effects 
of chemical fractionation and isotope-selective photodis- 
sociation, which would cause the 13 CO/ 12 CO abundance 
ratio to range between 1/70 to 1/35 (§3.5; van Dishoeck 
& Black 1988). 

The populations of the first eight rotational transitions 
of the 12 CO and 13 CO molecules were modelled. For each 
species the model was run for temperatures ranging from 4 
to 40 K, for densities ranging from 10 1 to 10 6 cm -3 , and for 
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column density per velocity ranging from N co /AV = 10 15 
to 10 18 cm- 2 /( kms" 1 ). Results for T k = 8, 22, 30, and 
40 K are displayed in Figure [n]. The observed range of 
values of the four ratios, Ri ± a{ (defined in Equations 
H through [ll]), are indicated by pairs of contours. The 
range in the (4-3)/(2-l) ratio is shaded. For reference 
the size along the line-of-sight, L, of CO emitting re- 
gions is indicated by straight solid lines. A linewidth of 
AV = 1 kms -1 has been assumed, whence the following 
relationship gives the size: 



L 

— 4.05x10 
pc 



N, 



(lO 3 cm- 3 ) 



1 -1 



8xl0- 5 

_ (17) 
Figure R2 shows the statistical analysis of the LVG re- 
sults. Using the method described, we have computed 
the probability distribution function, P(r)e?r, as a func- 
tion of N co /AV, n, and Tk- The greyscale images in 
Figure O each represent a slice through the probability 
density field made in the (N co /AV, n)-plane at a partic- 
ular temperature. The dashed contour is a slice through 
the the 3<r confidence surface if the (4-3)/(2-l) ratio is 
neglected and the solid contour represents a slice through 
the 3cr confidence surface when the new data are included. 
The crosses indicate the point in the (N co /AV, n)-plane 
where the probability density reaches its maximum value 
for the given temperature. 

The maximum probability density in each tempera- 
ture slice, Pmax(Tk), i.e., the value of the probability 
density at the location of the crosses, is printed on the 
plots. The probability field has two local maxima, one 
at T k « 8K [P max (T k ) = 1.92 x 10~ 3 ] and one at T k w 
20 K [P ma x(T k ) = 1.61 x 10 -3 ]. These correspond to the 
HDLT and LDHT solutions, respectively. The HDLT solu- 
tion is the most probable. This regime is very well-defined 
in the (N co /AV, n)-plane, with n — 10 4 5±a5 cm -3 , and 



N co /AV = 10 



16.6±0.3 , 



'( kms 



It is reassuring 



that we get the same value of N co /AV in the HDLT /LVG 
model as in the LTE model. The two models tell the same 
story: if the density is high then the level populations are 
thcrmalized in constant-opacity structures with T k « 8 K. 
The 12 CO line center opacities in the LVG model are: r = 
9.6, 14.4, 6.3, and 1.0 for the (1-0), (2-1), (3-2), and (4-3) 
transitions, respectively. In this model all 12 CO radiation 
is produced in gas which is at least marginally optically 
thick. 

As mentioned by FP the implied structure sizes in 
the HDLT case are very small, less than ~ 0.01 pc (~ 
2000 AU), which is consistent with the < 6500AU-sized 
structures seen by Pound et al. (1991) in high-latitude 
cloud MBM-12. This HDLT structure size is <;1% of the 
average projected CO size of HLCs (~ 1 pc; Ingalls 1999). 
Thus in this interpretation a typical HLC is comprised of 
an ensemble of ~ 10,000 HDLT CO-emitting cells (given 
an area filling fraction of unity) . 

Unfortunately, it is still not possible to rule out the 
LDHT solution completely, since the probability field has 
a second local maximum at T k w 20 K which is about 80% 
as probable as the 8 K maximum. We are, however, able 
to place a strict upper limit on the temperature in the 
CO-emitting regions due to the predicted size of LDHT 
structures. As mentioned in §1, the sizes of emitting struc- 
tures along the line-of-sight (LOS) may be much larger 



than their sizes in the plane of the sky. It is not plausi- 
ble, however, for LOS structure sizes to exceed the ~ lpc 
cloud sizes. For very high temperatures (T k > 30 K) the 
LOS size of emitting regions with nonzero probability in 
the model becomes unrealistically large. Unless most ob- 
served CO structures are filaments aligned along the line of 
sight, which is unlikely, temperatures greater than ~ 30 K 
must be completely ruled out. 

Again (see the LTE model) the 12 CO (4-3) data signifi- 
cantly decrease the size of the 3er surface in the parameter 
space. In other words, the range of possible solutions is 
decreased. Remarkably, including the higher transition 
data does not change the predicted physical conditions. 
We plot in Figure the maximum-probability values of 
P, (N co /AV), and n (the values at the positions of the 
crosses) . As for Figures lTJ through [IJJ, dashed curves rep- 
resent these functions when only the(3-2) and lower tran- 
sitions of CO are included, and solid curves represent the 
same functions when the new (4-3) data are used. Clearly 
the new data both limit the range of possible solutions in 
T k -space (Figure top) and also yield cloud conditions 
that are entirely consistent with those derived using the 
old data (Figure |l3|, middle and bottom) . 

5. DISCUSSION 

5.1. Physical Implications of HDLT Gas 

We have studied the J = (1-0) through (4-3) transitions 
of CO in translucent clouds. Motivated by the work of Fal- 
garone et al. (1991; 1992; 1998) and FP, instead of treating 
clouds as separate entities with distinct characteristics, we 
sought regularity in the CO emission. The linear relation- 
ships among the (1-0), (2-1), (3-2), and (4-3) transitions, 
as well as the constancy of the 60/^m/100^m IRAS surface 
brightness ratio, justify this perspective. The simplest in- 
terpretation of these results is that the emission arises in 
regions with the same physical properties, or "cells" , and 
that variations in the emission are caused by variations in 
the beam filling fraction of the cells. We call this the "cell 
hypothesis" . 

Two types of CO cell may exist: the low density, high 
temperature (LDHT) and the high density low tempera- 
ture (HDLT) cell. The LDHT cell is relatively space-filling, 
with line of sight sizes of order 50% to 100% of the cloud 
size. It is also in thermal pressure equilibrium with the 
interstellar environment. 

The existence of HDLT cells, on the other hand, is con- 
troversial. We have shown that cells with densities of or- 
der 10 4 5 cm -3 are the most probable source of CO emis- 
sion. This challenges the picture of HLCs described by van 
Dishocck et al. (1991), who argued against the possibility 
of gas with n ~ 10 4 cm -3 . The HDLT scenario also de- 
fies one of the primary assumptions of interstellar physics, 
which holds that clouds are in hydrostatic equilibrium. 
HDLT cells have a thermal pressure P/k ~ 2x 10 5 cm -3 K, 
which is an order of magnitude greater than the ambient 
thermal ISM pressure. Is this a serious problem? 

The interstellar medium is turbulent. Its structures are 
not the result of equilibrium processes. Recent simulations 
of magnetohydrodynamic (MHD) turbulence performed 
by Ballesteros-Paredes, Vasquez-Semadeni, & Scalo (1999) 
show that thermal pressure equilibrium is irrelevant, since 
in MHD turbulence the volume and surface terms of the 
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Fig. 13. — Maximum-probability values of P, (Nco/AV), and n in slices of constant (i.e., the values of the quantities at the crosses 
shown in Figure [13). Dashed curves represent these functions when only the (3-2) and lower transitions of CO are included, while solid 
curves include the new (4—3) data. The new data reduce the range of possible solutions, but do not alter significantly the most probable cloud 
conditions. Note that P m ax(T) is normalized so that the peak value equals unity. The absolute peak probability density value is 4 times 
higher for the current dataset than for the prior dataset. 



Eulcrian Virial Theorem are comparable. As a conse- 
quence, structures will nearly always deform under in- 
ertial motions, resulting in "clumping" on many scales. 
According to this model, in a turbulent cloud, clumps 
are continually forming due to density fluctuations and 
dissipating due to radiative cooling. Quasi-static config- 
urations are highly unlikely, unless compressed to sizes 
L < 0.01 pc, where self-gravity becomes important. Iron- 
ically, the HDLT cells are approximately this size. If the 
simulations of Ballesteros-Paredes et al. (1999) are correct, 
then the only stable configuration in translucent molecular 
clouds is the HDLT one. This important result needs to 
be confirmed. Is the CO emission produced only in sta- 
ble structures which are in LTE? This would explain the 
remarkable uniformity we see in the CO line ratios for a 
large sample of clouds. 

5.2. Limitations of the Method 

We have attempted in this paper to evaluate two com- 
peting models for the physical state of molecular gas in 



translucent clouds, the HDLT and LDHT models. Both 
van Dishoeck et al. (1991) and FP state that neither HDLT 
nor LDHT conditions uniquely describe the data (even 
though each group concludes by favoring a different con- 
figuration; see §1). We have assumed that clouds are 
macro-turbulent, and we have argued that the LVG ra- 
diative transfer model is an adequate approximation of 
macro-turbulence. By deriving a probability density field 
in the physical parameter space we have shown quantita- 
tively that HDLT gas is the more likely source of observed 
CO emission. 

Our analysis is limited, however, by the extent to which 
the LVG scenario is an accurate representation of the 
velocity field of translucent gas. The assumption that 
clouds are macro-turbulent is based on the extensive three- 
dimensional radiative transfer simulations done by Park & 
Hong (1995). Nevertheless, it can be argued that using the 
macro-turbulent LVG assumption has biased our results in 
favor of HDLT conditions, since tiny HDLT cells are more 
likely to satisfy the LVG criterion, because they are less 
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likely to shadow each other along the line of sight than 
the much larger (<~ 0.5 — lpc) LDHT cells. Furthermore, 
since they have larger values of N co /AF, LDHT cells 
have a much larger column density of radiatively coupled 
gas than HDLT cells with the same internal velocity dis- 
persion. Thus the macro-turbulent approximation should 
break down for the largest LDHT cells. A more complete 
study would vary the statistics of the velocity field between 
the micro-turbulent and macro-turbulent limits. This is 
beyond the scope of the current paper, and we leave it to 
future work. 

5.3. Testing for HDLT Cells 

If HDLT cells exist, then they can be detected directly in 
CO emission using a millimctcr-wavc interferometer. A re- 
solved HDLT cell would have a Rayleigh- Jeans brightness 
temperature in 12 CO (1-0) of T A ss 5K. An interferome- 
ter with a spatial resolution of 0.001 pc (for a cloud at a 
distance of 100 pc this corresponds to an angular resolu- 



tion of 2") would be able to detect such HDLT cells. If 
there is no detectable 12 CO (1-0) emission in HLCs at 2" 
scales, then the HDLT hypothesis is false. 
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